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IMPORTANT ANNOUNCEMENT 


The increase in the output of papers in the field of mathematical 
biophysics makes it difficult to insure prompt publication without an 
increase in the size of the journal. Therefore, the Bulletin of Mathe- 
matical Biophysics inaugurates the following service: 

Upon acceptance of a paper, the Editor, if necessary, will ask the 
author to shorten the paper to an extent dictated by the requirements 
of a reasonably prompt publication. The shortening should in no 
case reduce the paper to a mere abstract. Such a shortened paper will 
be published within six months or less. 

The unabbreviated original manuscript will be kept on file at 
the editorial office. Any person desiring to avail himself of the com- 
plete manuscript, may obtain promptly a microfilm copy of the latter, 
at the cost of 1¢ per page plus postage, by applying to the Editorial 
Office, 5822 Drexel Avenue, Chicago, Illinois. 

All papers in the Bulletin which have been thus shortened, will 
be marked at the end by the symbol MF, followed by a figure, indicat- 
ing the number of doublespaced typewritten pages of the unabbrevi- 
ated manuscript. 
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A MATHEMATICAL THEORY OF PROTOPLASMIC 
PROTRUSIONS. I 


HENRY E. STANTON 
THE UNIVERSITY OF CHICAGO 


The theory of cellular elongations is applied to a particular case of 
a spherical cell with a cylindrical protrusion under the conditions of equi- 
librium, and conditions for continued and arrested growth are examined. 


N. Rashevsky (1940) has developed and proposed a theory of 
cellular elongation derived from the mechanism of diffusion forces 
with primary reference to the phenomenon of cell division. We shall 
apply the same general theory in this paper to a particular type of 
cell configuration, wherein the main body of the cell will be consid- 
ered spherical and a cylindrical protrusion has been extended from 
this main body. In any metabolizing cell, the diffusing substance or 
substances will exert a drag on the cell structure, assuming that the 
diffusion takes place through some type of gel, and these diffusion 
forces together with osmotic pressure at the membrane and surface 
tension will cause deformation of the cellular shape under certain 
conditions. 

With these ideas, N. Rashevsky, and others have shown that the 
relative elongation in the z-direction is given by the formula, 


1d, _ 1 
Let P'3yV 
is 
aoa c{z cos(v,2) — 4[xcos(v,2) + ycos(y,y)]} dS (1) 
+ { [zZ, — 3 (aX, + y¥r)] as} 
8 
where 


(x,y,z) = coordinates of the surface elements dS, 


c= concentration at that point, 


X,, Y,, Z,=components along the coordinate axes of the forces 
~normal to the element dS, 
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y =a unit normal to the surface element, 
V = the total volume of the cell, 


u= the fraction of this total volume which is filled with 
semi-rigid structure, 


n = the coefficient of viscosity of this structure, 
R= the universal gas constant, 
T = the absolute temperature, 


L,= the length of the cell along the z-axis, along which 
the elongation is assumed proceeding. 


The components of the surface forces X,, Y,, Zy arise from the 
pressures, one of which is the osmotic pressure resulting from the 
change in concentration across the membrane, 


RT 
a een as ' 2, 
P=TF (¢ =) (2) 
where c’ is the concentration adjacent but outside, and the pressure 


arising from the surface tension, 
1 1 
SN Cee ee 3 
p y ( R, : ) (3) 


where FR, and R, are the principle radii of curvature at the point on 
the surface. These quantities are equivalent to hydrostatic pressures 
and are, therefore, directed along the normal to dS. 

The geometrical arrangement of the cell is shown in the diagram, 
and from this it is clear that the computations can be separated into 
three convenient parts, the main spherical body, the cylindrical sides 
of the protrusion and the end of the protrusion. We will consider first 
the evaluation of the integrals on the right side of (1) over the main 
spherical portion of the cell. 

It will be assumed that c. is a constant over the surface of the 
region, and that 6’ represents the angle subtended at the center of 
the sphere by the intersection curve of the protrusion, also considered 
a constant in these calculations. Then, for any point on the surface 
it is clear that 

z 


4p) 
cos(, 2) =, cos (vu) =< cos (vz) =. (4) 
4 


From (2) and (8), 
2y 


=e (6) —¢ ; ‘— _ 
D (C. 2) p ar 


M 
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which mean that the components of th : 
the sphere are e normal force at any point of 


FIGURE 1 
gat & 2 x 
Gs SS . — ¢@’ ees pee 
M (ce C's) =| r, erinrs (6) 


the other two components being derived by multiplying the bracketed 
expression with the proper direction cosine. Substituting these in 
(1) we have for the contribution of the spherical portion, the expres- 
sion, 


PaSRT eG 
=o f plese +) las 
g 2 


f RT ora ta U8 
= : FP (0) — | tae ty) as. 


This integral may be evaluated by using spherical coordinates 
and turns out to be 


; Slolie he 2y 
i,—2 1 cos 6 sin? ‘| (Gs — C3) +2] ; (8) 
Similarly by considering the cylindrical portion we have 
earners EE Cop — sed (Cop — C'2p) + =| . (9) 
2M M tT". 


At this point, it should be indicated that the assumption has been 
made that ¢.) remains constant over the entire length of the cylindri- 
cal surface. If we consider a cell like some of the yeasts which have 
very long protrusions, this would be a close approximation, especially 
if production of the substance still goes on in the protrusion. How- 
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ever, for short lengths as in pseudopods in amoebae, this approxima- 
tion may be doubtful. In this case, an expression which should be 
closer to the truth would be, 


C2 — Cip 
Cop — Co ra ae 


—2z (10) 
rs 
but this will not be considered at the present time. 

In finding the contribution to the relative elongation of the end 
of the protrusion, some approximations must be considered for con- 
venience. If the substance is considered as diffusing across the end in 
a direction always parallel to the axis, it makes no difference mathe- 
matically whether the diffusion forces and osmotic pressure are cal- 
culated over a flat end or over a hemispherical cap. The contribution 
from surface tension must be considered over a hemispherical cap 
mathematically. That the diffusion will not behave in this manner is 
rather obvious, since the substance will diffuse or flow normally across 
the boundary, but the approximation will be reasonably close under 
the simplifying assumptions made. If it is assumed that c,, is a con- 
stant over the spherical cap, the second method will be approximated. 

The calculations lead to 


BRD eT Spidey 
“om ° ae va (Cap = Cig) == = | (11) 


and we are now in a position to examine the relative elongation. 
Substituting (8) (9) and (11) for the expression in braces of 
(1), we have, after some simplification, 


le=anT's7 (71 +72 cos 0’) = 


Ich, vakT py 
rar aE = a0 |e" COS 6’ Sin?6’ + Copf"9?7", z Cg"? (15 
+ 7’, cos 6) 
2yM 


Seyi Ti [27’,? cos 6’ sin?6’ + 7’yr", — 2r”.(1', + 72 cos 6’) ] (12) 


2 , , , Ld 
ae [ (cz — ¢'2) 72° cos 6’ sin?6’ + (Cop — C'np) 1027", 


— (Crp — C'rp) 12? (7, + 7'2 COS 6') J 


which gives the relative rate of increase of length in the z direction. 
Of course, in this problem, there will be a contraction in the other 
two directions at right angles to this direction under the ae 
of cylindrical symmetry. 
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To simplify the expression let 
2My aRT u 2 
= BRT a : = 30M et 3a (13) 
and from the diagram, Figure 1, clearly, 
p= Po SiN:8 ss (14) 


The relative elongation then becomes 


Id; me 
Lo dt = Br’, sin ¢{—A[?r’, + 27’,(1 — sin 6’) cos 6] 
=i Ce sin Oe, FF Crp) To COs 6’ paar a (Ce ae C's oon, Cip te C'1p) (15) 


PTs COS OF (Crp 1615) 75 = OCdy “iC 55 = Cig t.Cas) Calin 


The terms involving a arise from the contribution of osmotic pres- 
sure, and if the enclosing membrane has infinite permeability (h) 
these will vanish. 

Elongation will take place when the right hand side of (15) is 
positive, and for equilibrium conditions it will be close to zero. Un- 
der the simplifying assumptions of infinite permeability, protrusion 
will occur when 


T° (C, —.Ci,)-Sin 6 cos, — Ar, 

+ [ (Gp. «,) oesine. — 2A (ii-<sinig') cos 6) 7 08 U6) 
Physical considerations of the mechanism of diffusion require 

that, for a cell which is producing a substance, 


ey ene San ye (17) 


and all quantities in parentheses involving the concentrations be 
positive... ; 

Hence (16) imposes immediately the restriction that r’, must be 
greater than the single positive zero of the quadratic in 7”, , and so 
determines a lower bound for the size of the cell in terms of 71, 0 
and the concentrations. It is more useful, however, to consider the 
restrictions imposed upon the other parameters, and so we arrange 
(16) in the form 

(7, + 7,:cos 6) [A — 72 (Cop — Crp) Sin 6] 
< 7’, cos 6 [A (2 sin 6’ — 1)-+ 1’, sin 6’ (Cy —-Cap) ] . > (12) 


We may consider separately the following possibilities: 
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Lee At rs (Co Cap) SI 
>A(2sin 6 —1) +72 sin 6'(C. — Cop) > 0; 
LIS A 15 (CG) Sa 
> A(2sin 6’ — 1) +7. sin 6’ (C2. — Cop) 5 
III. A(2sin 6 —1)+ 72 sin 6’ (Co — Cop) 
> Ara (Cap 1) Sin oe 0; 
lV. 0>ACsin 6) 1) eee sie (Cc) (1 
> Ai To (Gp Cis) SIe G 
Vi Ait (Coe Cap RSI 
>0>A(2 sin 0-1) +72 8in O(c. — G5)G 
VI.“ A(2. sin. 6 — 1) 4+ rs sin? (cs) 
oS OSS Ate (Can Cag) SLL 


Cases I and V cannot satisfy (18) for any positive value for 7’; 
and consequently represent conditions of insufficient production and 
no consequent growth. Cases IV and VI will meet the requirements 
for any value of 7’; , and case II also, provided 7’, is sufficiently large. 
These three cases represent necessary conditions for continued 
growth. In contrast to the abundance of these, case III is the only 
one which establishes an upper limit for the length, 7’, . 

By means of (18) and (19) various relations can be established 
among the parameters. We suppose the cell to be producing the dif- 
fusing substance so that (17) holds. All of the conditions, II, IV, VI 
and III, require 


Co > Cip (20) 
automatically, but if in addition ec > ¢., they also demand that 
sin é@, <4 (21) 
or from (14) 
eo. ae (22) 


In case II under these conditions, it follows that 
Ts Cos 0 2A (l= sin 6)))— 7, sim te, — C35) } 
, ° , < La < (oe) 
Tz SIN 6’ (Cop — Crp) — A 
= sin ¢ S$ ——_—____ (23) 


A 2A 
1'2 (Cop — Crp) 2A aT, (C3~"C35) 
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Cy tl 183 23816: ie os eae Cae 

2(1 — sin 6’) ” 2a 
Case IV similarly gives 

OQ Ser <5 
Sai 2A 
sin 6’ > DAC (Soe: (24) 

Gor ig (1-02 Sin), 

2(1 — sin 6’) 
The middle inequality together with (21) requires 

A 

Cea Ci, 

The only difference between cases IV and VI involves the con- 
centration ¢,, which is limited in range such that 
A(2 sin 6’ — 1) 


Cop Cie 
r’, sin 6 ftir 


D 


2p 


Te. 


(25) 


rain o as 
The final case to be considered is that permitting only limited 
growth, case III, which gives the relations 


Te COS 2A (Sin 04-41) 775 Sin 6 (6, — Can) | 
A- we sin @ (Cs, aa Crp) 
2A 


sin 6° > ——______—__— 
2Ae T's Co Bia Crp) 20) 


Cael —-ZSitb ) Cus 
i 
Zl. ==)'Sin.6.) 

These relations are necessary conditions for growth to occur. 
They all seem to be very similar in the terms involved and the various 
dependencies. The concentration c., on the sides of the protrusion 
will be directly dependent on the rate of production in the protrusion, 
and it will increase as the production increases, particularly when the 
protrusion is very long, or relatively narrow so that diffusion from 
the main cell body becomes relatively unimportant in the protrusion 
except in the neighborhood of the base. For this reason, the relations 
involving ¢:) may give some indication as to the extent of activity in 
the protrusion. 

With this in mind, one would expect that when only limited 
growth is permitted, the activity in the protrusion would be a mini- 
mum, whereas higher activities would cause the growth to continue 


0<ri< 
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indefinitely. Of the cases allowing this, II permits the narrowest pro- 
trusion, but because of the lower limit for 7, , the parameters must 
be adjusted so that one of the other cases applies first, whereupon 
case II may be followed. It might seem apparent that with sufficiently 
large cells which are following this pattern, a point might be reached 
where the protrusion had vanishing diameter. This is, however, not 
the case, for by considering (14) and large 1’, we find 
2A ae Co + Cip 


” — 
Lae ’ Cop — 


28 
ty) 5 (28) 


which gives a value of the radius of the protrusion depending only on 
the surface tension constant and the concentrations. 

In the calculations, no assumption has been made of constant 
volume, and this picture requires that fluid diffuse into the cell suf- 
ficiently fast to: permit growth, and the phenomena are, therefore, 
only long time ones to which this analysis would apply. This approach 
gives considerable simplification in the use and interpretation of Bet- 
ti’s formula, since 7’, can be considered constant. This means that L, 
will be a linear function of 7’; only and the relative elongation rate 
can be thought of as involving only this quantity as changing with 
time. Also lateral contractions which are implicit with the assump- 
tion of constant volume do not enter. 

If constant volume is assumed, which includes the more general 
case, L, would be a function of 7’, and 7’, and from physical considera- 
tions, only limited growth would be possible since the main cell body 
would decrease in size aS growth progressed. Also lateral contrac- 
tions should properly be calculated. This extended problem might 
well apply to numerous cells and growths which cannot be treated in 
this simplified discussion, such as amoeboid movements. The present 
paper may have applications in yeast and nerve growths. 

The writer wishes to thank Dr. A. S. Householder for valuable 
discussion concerning the mathematics and conclusions, and Dr. N. 
Rashevsky for suggesting the problem together with ideas about its 
application to the theory of growth. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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A THEORY OF ELECTRICAL POLARITY IN CELLS: II 


ROBERT R. WILLIAMSON 
THE UNIVERSITY OF CHICAGO 


Orders of magnitude and relations of the constants in the equations 
of the polar cell are discussed. Some considerations on the mechanism 
of polarity are presented, and a more general mechanism discussed. It is 
Sen that polarity can exist without catalysis, either MEER or in- 

ibiting. 


In a preceding paper (Williamson and Bloch, 1942) equations 
were derived determining the polarity of a cell maintained by diffu- 
sion and electrical forces. The diffusion forces acted in two ways, i.e. ; 
(1) they acted directly on the catalyst particles to produce a polarity 
if the particles inhibited the reaction responsible for the diffusion 
forces, or to prevent a polarity if the particles promoted the reaction, 
and (2) they separated ions of opposite charge produced by the re- 
action by a differential in diffusion rates. The second action leads to 
an electrical polarity of the cell which influences the distribution of ' 
the colloidal particles if the catalyst particle has a:charge. In the case 
of an inhibitory catalyst the diffusion forces of the first kind always’ 
tend to establish a polarity, while in the case of a promoting catalyst 
they will oppose a polarity. The electrical forces. will either aid. or: 
oppose a polarity depending on the signs of (D, — D-_) and » (the. 
charge born by the catalyst.) . a 

The equations determining polarity are a lipaieon and Bloch, 


1942) 


- : . 0 aN 
ploa by ) 
Ge an, ERE patos oustnnn tae 
Q.=Q. — an, ies St 
3Nv - (2) 
a= DE 
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tea pal 1 ut Lied 
coeds Daa) ala a) 
_ 56 a7? N » & 
B= K - 


If we let 7 = 10° cm, 7, = 10° cm to 10° cm, D,— 10° cm? sec’; 
D_= 10-5 cm? sec”, K = 80 we obtain as the upper and lower limits 
of M 


M < 410"a + 5.10"a — 1.7-10+ 
< 5-10*%a — 1.7-10-, 

M > 4-10?a — 5-1078a — 1.7-10°° 
> 5-10a. 


Thus we see that M may cover a very wide range. 
If |My| is large 


My My 
é 2 —e 2 
di ana ie 
AS a 
i.e., either 
nm, = 2nandn,=0whenM>0, 
or 


N, = 2nandn,=0whenM <0. 


When My is small we can expand tanh My/2, and take only the first 
two terms. This gives a better approximation for |My| 2 1 than the 
first three terms, and a good approximation for |My| < 2. 

The equation then has the stable roots 


_ 2 /8nM—3\ 
¥~"“M\ nM ). 


and the unstable root y = 0. When a > 0 it is necessary that 
Qo = an, Q, 2 an, 
to prevent negative reactions. Then 


eg ee 


) an. 


The non-null real roots will exist only if nM > 1. This is equivalent 
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to the condition 


d 

— tanh bd) >—, 

dy 2 |, 2n 

@ 

which is geometrically necessary, as seen by examination of the 
graphs of tanh My/2. 
____If we neglect fy/RT , M = aB, and the upper and lower limits 
indicate 


B& 
— 5-108 S B= yt + — = 5.10% 
vies iT = 5-104, 
If we make My < 2 so that we can use the expansion, we have 


1 a(n, — Ne) 
= >= > ———_. 
Q. > an 3 > 5 
Thus if Q. ~ 107 moles cm sec", then an < 10-7 moles cm~ sec-. 
Furthermore, since the volume of the cell is of the order of 4-10-° em?, 
and as there must be of the order of 10 to 1000 particles in the cell if 
the distribution is to have any statistical uniformity, we have 


10 ne 10° 
410° 4-10’ 


or n & 2.5-10° to 2.5-10"' particles em at least. This corresponds, if 
each particle bears one molecule of catalyst, to a solution of the cata- 
lyst of 410-2 to 4-10-° moles liter-*. This is at the lower limit of dilu- 
tion of the most potent biologicals, and therefore we would expect 
either a higher » or more catalyst molecules per particle. There is 
another factor, however, which alters these considerations. The vol- 
ume of the cell is 477,°/3 , while the volume of each particle is 427,3/3 . 
The relative volume of active protoplasm left is then 1 — 4znr,3/3. 
We must assume then that the catalyst concentration is, in terms of 


active protoplasm 


NCW 


ett Oe 

—N 1— 4nnr,2/3 
However, because the catalyst molecule is screened by the inert mass 
of the particle on one side, the effective concentration is roughly only 
one-half of this since only half as many molecules can get to the cata- 
lyst in a given time interval as when the catalyst is freely dissolved 


in solution. ' 
If we revise the derivation consistently with the assumption that 


Ce moles/liter. 


104! ELECTRICAL POLARITY: f'N’ CELLS 


the particles are impenetrable to the metabolites, we must -intro-- 
duce two new changes in the original derivation. The diffusion co- 
efficients of the metabolites. will, be altered. by the presence of these 
particles. Einstein 1906) has: shown that a suspension of rigid 
spheres of, uniform, size, whose common. radius is greatly exceeded by., 
the mean distance between the particles and the smallest dimension, 
of. the shearing liquid, has a viscosity, 


ni = No (1 + 2.5 vn), 


where vn; is the fraction of the volume occupied by the dispersed par- 
ticles and mp is the’ viscosity of the fluid alone. Also Smoluchowski 
(1916) has shown that if the particles carry an electric charge the 
viscosity ix further increased by an electro-viscous effect thus 


an ‘ 1/K¢ 
Peto eapenll 


Pee res ee age hd = specific conductivity, | 
i ‘K = dielectric constant, © 
¢ = electrokinetic potential. 
If we let 


a= 25/14 mace 
ed ANo 221, 


ni =no(1 + avni). 


Einstein’s expression for the diffusion coefficient is 


| kT 
sag ees AiG te 6 ant 
If we. tates the the vires and expand we have 
a OE at if kT ; kT’: ; 
Gay(ltavni)r 6anor (1 —'o Bite) 
=D Alarm). 


Because of the many crude assumptions the above can be peste 
only a first approximation. 

The second change required in the original derivation’is the fol- 
lowing. ‘Tf Q> ‘is the production rate of a metabolite in protoplasm 
without granules, we have, for protoplasm containing granules, 


Wold aeadtpeniene of! le Qy =e Quins Vane yo t 


ie ue 
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and if the particles bear a catalyst, 
Qi = Qo(1— uni) + bni = Qo — (Quy — b) nj 
= Qo — an; , 
where 
a—Q,v —)d. 


This equation is identical in form with the previous equation, 
but has a more general physical significance. It now becomes appar- 
ent that a polarity can occur even without the catalyst. We should 
expect self regulating polarities resulting from this mechanism in 
cells having no structure more specific than some granules, or a 
vacuole. 

The solution of the equations resulting becomes extremely in- 
volved because of the non-linearities introduced. They must reduce 
to the equations of Williamson and Bloch (1942) as a first approxi- 
mation, however, and for the preliminary applications to biological 
problems the older equations will do very well. 

This work has been aided in part by a grant from the Dr. Wal- 
lace C. and Clara A. Abbott Memorial Fund of the University of 
Chicago. 
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GREEN’S FUNCTIONS IN BIOLOGICAL 
POTENTIAL PROBLEMS 


ALVIN M. WEINBERG 
THE UNIVERSITY OF CHICAGO 


Electrical potential problems encountered in biology differ from 
those usually considered in electrical theory first, because the membranes 
of tissues satisfy a non-linear relation between current flow and polari- 
zation, and second, because the interior of the tissues are not equipoten- 
tials. A Green’s function suitable for discussing such problems is de- 
fined, and a cylindrical illustration of such a function is discussed. 


Potential problems arising in the electrophysiology of nerve and 
muscle are of a type not usually encountered in electrical theory for 
two reasons: first, because of its finite internal conductivity, the in- 
terior of an active or injured fiber is not strictly equipotential; sec- 
ond, because the membrane resistance is generally a function of the 
membrane current flow, (Cole and Curtis, 1941) the boundary condi- 
tions are non-linear. In the present note a Green’s function suitable 
for the solution of such problems will be defined and a cylindrical il- 
lustration due to H. Weber will be discussed. The results are at pres- 
ent primarily of mathematical interest, and are rather closely related 
to some results first stated by H. Helmholtz (1853). It is hoped that 
they may constitute a ground work for future attempts to analyze 
bioelectric potentials by the methods of potential theory. 

1. The Green’s function. 

Suppose u, and u; are the steady state potentials outside and in- 

side a tissue covered by a closed polarized membrane, S (Figure 1). 


P (xy,2) 


FIGURE 1 
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Then ue and u; are harmonic functions which satisfy, in general, a 
non-linear boundary condition at the membrane: 


Oo eee eee ; 1 
a k; A F(u. — uw) on Sy (1) 


2 [Ole OU; 
um = a(S) =o. (5): 

where k, and k; are the conductivities of the external and internal 
media, 2, and n; are the outwardly and inwardly drawn normals on 
S,and F(u. — u;) is some function of the membrane potential dif- 
ference which determines the membrane current flow. The exact 
form of F' has been determined experimentally by K. S. Cole and H. 
J. Curtis (1941) in the squid giant axon. Prior to this work, most 
writers on the theory of the propagated impulse had assumed a piece- 
wise linear relation between the membrane current and the membrane 
potential difference of the general form 


é 


or 


F (tte — us) = TE (te = Us — B) (2) 


where FE, the “membrane E.M.F.”, is constant except in active re- 
gions, where it is zero (Rushton, 1937). The constant h is the ratio 
of the external conductivity to the membrane conductivity per square 
centimeter; it is a constant only for very small membrane current 
density. To complete the specification of the problem, we shall assume 
either that the vessel, V , in which the tissue is immersed is an insu- 
lator, in which case 


OUe 
ONe 


=0 on VJ, | : (3) 


or that it is a “diffuse ground”, in which case 
Ue=0 onV. (4) 


Both of these experimental situations arise in practice. 
Consider an exterior point P(x,y,z) from which the distance Pe 
to a variable point P’(é,7,¢) is 
Park Odea) lane, oe) la mel (8 ed (5) 


The reciprocal distance 1/p, has a pole in the exterior at the point P, 
while it is a continuous function everywhere in the interior. Conse- 
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quently from Green’s theorem 


1 Oe 0(1/pe) 
ae Ue ‘ |eo 
pe ON; On; 
am Me OUe 0(1/pe) ¥ 
pe ONe ie ON. | Oe 
=e aus pale, 
ey Pe ea ONe iets ) 


where w and w’ are surface elements on S and V respectively. 
Let us now introduce two harmonic functions, H, defined in the 
exterior, and H; defined in the interior, which have the properties: 


OG? k; 0Gi? 
ON Pee Fic On; 


(6) 


and 


—=(G.°> Gi") aon 55 (8) 


OG AN. OU, 
—=— if 
On; Vv ONe 
where v is the surface area of the insulating vessel (cf. Kellog, 1929, 
p. 246) and 
G76 23254, 0))= Wet HE, 4,2); 
(%,Y.52) 


=0 onV, (9) 


Gf=01f ue=0; or 


| in exterior, (10) 
(€,7,0) 
and 
ee (lee eh, = 1 pact tis (ph), 
(«,y,2) in exterior, 


(11) 
(€,,¢) in interior. 
Then, since H, and H; are both harmonic, 
1 f f OUe oH, 3 
== x ‘ H. an; i the an, @ 
(12) 


il f\ OUe oH a 
mu oe, ow, 
oH; 
sett ier Jaw. (13) 
ana Oe nt, 


and 
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If we add equations (6), (7), (12) and (13), and use the conditions 
(8), (10) and (11), we find 


nawaa=af Sle + EO a, 


=, (Gi° 5) G.°) (Ue — Ui) | dw (14) 


fff now 
Th E 5 ie ae Oe Ow 


The last integral vanishes if the diffuse ground is used; otherwise 
(insulating vessel) its value is 


1 ; 
to => | [ (En do's (15) 


this represents the average value of the potential over the containing 
vessel. 

The first integral may be simplified by substituting equations (1) 
and (8) and using the fact that d/on; = — 0/on,. The final result is, 
for a diffuse ground boundary condition 


wen. =75 J Ion * “Gr 


M [Gots Vek So Ya) Gye Ye eee yy ere, oy 
and for an insulating vessel boundary condition 
eS, 
te. Ys 2) = ies 0, (= “YJte. (@,y,2;&,,C) 
(17) 


me OFAN Cel er tet aney ee ae Meee ree. 


The additive constant a, in the insulating vessel case is required 
since only normal derivatives and differences between u, and u;, 
neither of which can determine a zero of potential, enter in the speci- 
fication of the problem (cf. Kellog, 1929, p. 246). In the diffuse 
ground case (t%, = 0), the zero of potential is determined uniquely 
by the fact that the vessel is assumed to be a perfect conductor at 
ground potential. 


Expressions analogous to equations (16) and (17) can be given 
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for the interior potential; for example, in the insulating vessel case, 


_ ke eh 02; of CU; pee ie . 
ACS yee) aie Ee Zeal i ea ® ee 
s 


where the G;' and G.’ are interior Green’s functions defined over the 
following ranges: 


: (~,y,2) 
Gee pesc 4, C) in interior, 
(25a) 
. (oS 5 2) in interior, 
GE eon yc) 
(Exc) in exterior, 


and satisfying equations (8) and (9). In biological applications, the 
interior problem is usually of less interest than the exterior problem. 
This is because the volume of many irritable tissues (e.g. individual 
nerve fibers) is so small that accurate interior measurements are very 
difficult (see, however, Graham, Carlson, and Gerard, 1942). 

As an example of the sort of biological conclusion which can be 
drawn from equations (16) and (17), we may prove the rather self- 
evident statement: the external potential of a resting tissue is every- 
where constant ( = ty). For, by definition, in a resting tissue the 
membrane current, 62,./en; , is identically zero. Consequently, (0) 
is identically equal to LZ, the resting membrane E.M.F. (a constant), 
over the tissue surface, and, from equations (17) and (18) 


—k,E 0G ie 
Ue (%,.Y, 2) = @ + thes 
wes 


On; 
= Uev 
since G;° is harmonic in the interior and the surface integral over a 


closed surface of the normal derivative of a harmonic function is zero 


by Gauss’s theorem. 

It may be mentioned that practically all of the discussion of the 
“dipole” hypothesis (Wilson, McCleod and Barker, 1933) can be put 
on a rigorous mathematical basis by application of the formulae of 


this section. 


2. Weber’s two-phase kernleiter.' 


1 Weber’s one-phase kernleiter, in which the interior is a perfect conductor, 
has been considered before (Weinberg, 1941). The sign of hf in the discussion of 
the infinite external medium should be reversed in that paper. 
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An example of a Green’s function of the type discussed in 1 was 
given by H. Weber (1873) in discussing the stationary current flow 
around an infinite cylindrical kernleiter (radius b) surrounded by a 
concentric glass cylinder of radius a(> b), the intervening medium 
being filled with an electrolyte. The interior and exterior Green’s 
functions-are defined by 


aIGee" A oG,°* o°G,°? 02G.°* 
or’? Te eote Op? oz? 


=0 (OIC Te <= @); 


oPGiSt Ah 0G,°* il 02G,°* i 0G 4%" 
or’ yr’ or’ ye Op? 022 


2 (0<r' <b), 


with the following boundary conditions: 


0G,°* 
Gia ——(7.4" = () at 2’ = 0; a {) EX ae hee os 
or 
oC Ringe kj CS aren 
or Ty er ee — Gi® atr’=a 
0G. 
permet re earn CAG a atz=2z', 
0G;°* 
carck ce, eae atz=2', 


where 7’, ¢’, 2’ are cylindrical coordinates of the field point, f bs OL 
the source point, superscripts refer to the position of the source point, 
6 is the Dirac delta function, and h is the ratio of external to mem- 
brane conductivity per square centimeter. The solution of this prob- 
lem as given by H. Weber is 


Dizee went, 
Gilt, $5%575h,2 ) 
b<r<a 
1 
Teena pt ae Ia inl) 


X cos n(d — ¢') e-8l#-*"1 
OS 0 


Gy Fe, 22 poe ) 
b<mr<a 
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1 n=00 En 
SRE Dana Oe (r10 


X cos 0(4 — ¢’) e-9!7-*'1 


and 
Bar ia 
Galt Gees 40 5 2.) 
Oi 7 <0 
1 N=00 En 
2 oN, 00 9) fn (7 10.) 


a hee kes e 


x cos n(¢ tz} ¢') e-912-*"| 


Om gy -<ah 
Gono eit, $572) 
O< 1D 
aoe! 1 co En , 
ge k;? n=0 =on,, on” medal!) 
X cos n(¢ — ¢’) e-8!?-*'1 
where 
fn(1,0) =Jn(O7) Yn (00) — Y,(07) Jn (0G), 
dn(’ , 6) =adn(67), 


a = Ke[In(00) Yn (00) — ¥n(0b)In' (6.4) ]/[hkidn' (6b) + Kedn (6b) ], 


1 ™b 1 a 
Nw =— | TP (OF) Oi if rfn2(07r)adr, 
Cad) ev, 


f= oe 2, 


&—1, 


and the characteristic values 6 are roots of — 
hok:=Ki[In(9 6) ¥'n(04) — Yn(0 db) In (64) | /[I'n (6 6) Yn (6 @) 


— Y',(66) J’n(8@)] — ke In(6 6) /J'n (8 6). 
The symbols J, and Y, stand for Bessel functions of the first and sec- 


ond kinds respectively ; primes represent derivatives. 


The characteristic equation reduces to 
h6= [F,(60 0) Y',(6@) — Y,(6:0) J'n(90)]/ [F'n (00) Y'n(@ a) 


YR (OD) F'n (Ga) ] 


when k; > > ke; this ‘case was discussed in a previous paper. On the 
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other hand, when k, > > k;, 


hoe s'.(8 b) +J,(0b) =0° 
this is the characteristic equation for an interior cylindrical problem, 
and its roots have been tabulated (Carslaw, p. 321). 

An application of this Green’s function which is of some biologi- 
cal interest is the following: consider an infinitely long kernleiter 
whose surface carries a uniform polarization, FE , from z=0toz= 0, 
and no polarization from z= — » toz=0. Such a partially depo- 
larized kernleiter may be taken to represent a nerve fiber in which an 
injury current is flowing. The boundary conditions at r = a are of 
the form 


Ue OU; ie ake Oa 
oer k; an; h (Ue Ui EF) rae Ves 
Ke 
= 5, (Ue — Ui) Ze 


To determine the external potential we may substitute these con- 
ditions into equation (14) ; we obtain 


E (oa) QT 
Ue(1, d, 2) = | ff G(T. do 9255 <2) 


—Gi (1, b 23 ed ae) | 0 des 


In the integration all Bessel functions of order greater than zero drop 
out: the result is 


aH 1 fo(@, 6) ho (a, 6) 
Aes = — YF ~~ —q- — um rake 
Ue(1, 2) 3 ae yar |roo, oe z2<0, 
ak il Ee 0) do(a, 6) 

aE are Pee EON Tl OS, 
For sufficiently large values of z only the first term in the expansion 


need be kept; hence, to good approximation, ignoring, however, ques- 
tions of convergence of the infinite series (cf. Jaeger, 1941) 


Ue (1 ,Z) =const X fo(7 , 6) e-% Z0) (18) 
Ue(7 , Zz) = const X fo(7, 6) (1 — e-%) z>0; (19) 


in like manner the interior potential is found to vary exponentially 
for large z; consequently the membrane potential difference will also 
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be an exponential function of the distance from the injured region, 
the length constant being given by 1/6, . 


An experimental study of the radial variation of the potential 
and current density around a passive iron wire conducting an impulse 
has been begun by Prof. F. L. Verwiebe and Mr. J. Lorenz. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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FURTHER CONTRIBUTIONS TO THE MATHEMATICAL 
BIOPHYSICS OF VISUAL AESTHETICS 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


In previous papers cases were considered in which a visual pattern 
consists of a relatively small number of relatively strongly excited ele- 
ments. Those cases are of interest in the theory of visual perception and 
aesthetic rating of artificial man made patterns. By an extension of ‘the 
theory so as to take into account the finite thresholds of the inhibiting 
fibers a theory of aesthetic ratings of pattern consisting of a very large 
number of elements is outlined. This type of theory is applicable to nat- 
ural patterns, such as landscapes. 


In preceding publications (Rashevsky, 1938, 1940), we discussed 
the total central excitation resulting from the contemplation of some 
simple geometric pattern and by connecting that total excitation with 
the sensation of pleasure it was possible in some cases to predict theo- 
retically the experimental rank order aesthetic ratings of some pat- 
terns. The excitation of each element of the pattern was considered 
as sufficiently strong, so that the thresholds of the inhibiting cross 
connecting fibers could be neglected. Under those conditions an opti- 
mum total excitation EH is obtained for a definite number of elements. 
A pattern containing too many elements has a small E' and therefore 
a small aesthetic rating, unless certain conditions of symmetry or of 
periodic repetition decrease the “effective” number of the elements. 

Observation and experiment tend in general to confirm these 
conclusions, so far as geometric artificial or so to say “man made” 
patterns are concerned (Rashevsky, 1940). But when we consider 
the case of a beautiful scenery or a landscape painting, the simple 
theory is clearly inapplicable. For simple geometrical patterns we 
considered as elements such things as straight lines, angles, lengths, 
ete. Neurobiophysical structures were outlined which result in the 
perception of the above mentioned things as “elements”. Similarly 
we can conceive of neurobiophysical mechanisms, which will corre- 
spond to the perception of a given color or shade as an “element” of 
the pattern. But a scenery or landscape painting consists of a tre- 
mendous number of such elements. The simple theory, applicable to 
elementary geometrical figures would predict complete mutual in- 
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hibition of such a large number of elements, which therefore could 
not be expected to have any aesthetic value at ail. 

Things are different however, if we consider the finite thresholds 
of the inhibitory fibers, leaving the rest of the neurobiophysical 
scheme unchanged. Equations for that case have been actually devel- 
oped (Rashevsky, 1938, chap. xvii), in another connection. 

Considering as before a set of n parallel excitatory neuronic 
chains corresponding to n peripheral stimuli or elements cross-con- 
nected by inhibiting chains and using the same notations as in loc. 
cit., we find that if 


hs 
LAS SALI rrr (1) 
(equation 56 of loc. cit.) then the inhibiting fibers remain unexcited 
and a central excitation results no matter how large the number vn of 
elements of intensity S,. As a matter of fact the total central excita- 
tion increases under these conditions with increasing 7”. 
If, however, 


i hs 

See ae (2) 

(equation 59 of loc. cit.) then the inhibiting fibers come into play and 
this results in a complete central inhibition when 

P OQ I 1 S Biel ae hy 
gti Gti IRL 4 St 
Only when h; = 0, inequality (2) is satisfied so soon as the first in- 
equality (1) holds. 

We may thus perceive a pattern consisting of a very large num- 
ber of elements and derive an aesthetic satisfaction from it as long 
as the intensity of excitation of the different elements is not too 
great. The equations derived in loc. cit. hold for the case when all ” 
elements are excited with equal intensity. But this can be readily 
generalized for different intensities of excitation, S;, provided the 
S;’s satisfy inequalities (1). Under these conditions a reduction of ” 
results in a reduction of EF and hence of the aesthetic value. On the 
other hand, keeping » constant but increasing some of the S; suffi- 
ciently strongly, so as to satisfy inequality (3) may produce two ef- 
fects. If the number of the strong S,’s is small and the strength of 
these S;’s is not too large, then although they will to some extent in- 
hibit the remaining weaker elements, yet due to their own increased 
intensity of excitation, the total intensity E may even increase. If 
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however many S;’s are strongly excited, this results in a complete in- 
hibition. 

The situation could be studied by using photographic pictures of 
the same landscape, showing different amount of details. Common 
experience tells us that an underexposed picture, lacking in details, is 
not so pleasant as a properly exposed one. An overexposed picture 
may both lack in details, as well as have too strong contrasts, and 
will not be pleasant either. 

The difficulty for a quantitative experimental study lies not so 
much in the production of pictures with different number of elements 
as in the actual counting of these elements. The following suggestion 
may. perhaps be used. 

When a photograph of a landscape is scanned on a facsimile 
transmitter, the variations of shade are transformed into electrical 
impulses. By taking a continuous record of the current, which leaves 
the transmitter while the picture is being scanned, we shall obtain a 
line which has the more ripples or peaks, the greater number of ele- 
ments the picture has. The number of individual ripples on the line 
of the recorded current may be used as a measure of the total number 
of elements on the picture. The record of the current produced by 
scanning a perfectly uniform sheet, with no elements, will give a 
horizontal line. 

The same method could be used in principle for the production 
of pictures with different number of elements, starting with a “stand- 
ard” picture. At the receiving end of a facsimile transmission out- 
fit the current impulses are again transformed into photographic im- 
pressions of different shades of black and white. If, by suitable con- 
densers or inductances: we distort the current curve in such a way 
that the smaller ripples on the curve disappear and feed such a dis- 
torted current into the receiver, we shall obtain a picture which has 
fewer elements. On the contrary, if by proper means, we distort the 
current curve in such a way, as not only to preserve all the ripples but 
to accentuate them, we shall have the case of increased S;’s. The pic- 
ture thus obtained will have the same number of elements as the 
“standard” picture, but they will be much more emphasized. 

By such a method one could obtain a set of pictures with different 
numbers of elements, and a corresponding set of current-records, 
which will give us directly the number of elements. The aesthetic 
values of the pictures could then be determined by standard psycho- 


physical methods. ; 
Individuals with very high hf, and A, will enjoy much stronger 


1 This has been suggested to the Author by Mr. H. E. Stanton. 
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S;’s , in other words much more pronounced contrasts in either shade 
or color. On the other hand individuals with very low h,; will either 
prefer very simple patterns or like complex patterns with extremely 
vague contrasts. Thus we have a neurobiophysical interpretation of 
differences in tastes. Other neurological parameters also enter of 
course into the question as seen from equations (1), (2) and (8). 
Some peculiar characteristics of modernistic paintings with either 
rather “crying” colors or extreme ‘‘vagueness” may perhaps be inter- 
preted in the above way in neurological terms. 

Similar considerations may hold about sound patterns. High 
thresholds of inhibiting fibers would result in the preference of 
strong, harsh and sudden sound combinations. 
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SOME OBSERVATIONS ON THE SIMPLE NEURON CIRCUIT 


WALTER PITTS 
THE UNIVERSITY OF CHICAGO 


A new point of view in the theory of neuron networks is here adum- 
brated in its relation to the simple circuit: it is shown how these methods 
enable us to extend considerably and unify previous results for this case 
in a much simpler way. 


I. Introduction and Definitions 

The sequent remarks are intended primarily as a prolegomenon 
and an illustration of a more general theory of neuron networks to 
be later redacted. Herein, as in the more general discussion, we shall 
adopt the linear model of excitation as employed by A. S. Householder 
in four previous papers in this Bulletin; (Householder 1941a, b, c, 
1942) our point of view will be somewhat different, however, in that 
we shall be primarily concerned with non-steady-state activity and 
the conditions under which a steady-state may be attained, rather 
than in the ascertainment of equilibrium points per se. We shall 
nevertheless find that this point of view has interesting consequences 
for problems of network statics. 

We may commence by defining the total conduction time of a 
fiber as the sum of its conduction time and the synaptic delay at the 
postliminary synapse: we shall suppose, as can be done without sub- 
stantive loss of generality, that all the total conduction times of fibers 
of the circuit GC in question are equal: and we shall measure time so 
that this quantity is unity. We shall consider explicitly only the case 
where all fibers of GC have remained unstimulated up to t = 0, a con- 
stant stimulation 4; + h; is applied to each synapse s; in the interval 
t «(0, 1), and the external stimulation assumes thereafter a steady 
value of s;. We shall indicate desinently how all other cases are sub- 
sumable with inessential novations under this one. . 

We shall employ the same notation as in (Householder, 1941a, b, 
ce, 1942) with ‘A’ for the product of the activity parameters of C and 
the abbreviation y; = o:\”), where there are n fibers in C. In addition 
to these, we shall find the following notations of value. 

The excitation-pattern of C may be described in a matrix E, of 
nm rows and an infinite number of columns, each of whose elements 
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€,, represents the excitation at the synapse s, during the interval 
(s,s + 1). The successive entries in the excitation matrix E' may be 
computed recursively from those in its first column—these are the 
quantities 4,—by the following rule, whose validity is evident: Given 
the elements of the p-th column, compute those of the » + 1-st thus: 
if the element ¢;, is negative or zero, place ois, in the 7 + 1-st row and 
p + 1-st column, or in the first row of the » + 1-st column if i= 7. 
Otherwise put oj + @€i, in this place. We shall say that C is in a 
steady-state during a series of n intervals (s,s +1),---,(s+u-—1, 
s +n) if, for every p between s and s + 1, the p-th and p + n-th col- 
umns of E are identical. If s is the smallest integer for which this is 
the case, we shall say that the steady state begins at the interval 
(ee SL). 

It will be seen that the construction of the matrix H implies that 
its infinite diagonals—where we take a diagonal to start again at the 
top of the succedent column whenever it reaches the last row of E— 
are wholly independent of one another, so that if we know the start- 
ing point of a diagonal of #,, we can calculate the entries along it 
uncognizant of any other values in the matrix. Physically, this of 
course means that the activity in C can be regarded as composed of 
wholly independent impulses, commencing originally at a synapse §$; 
with a value 4;, and journeying around C in irrelation to the im- 
pulses beginning at other synapses. We shall find it convenient to 
adopt this standpoint, and consider only the case of a single impulse, 
so that the complete solution must be derived by combining the re- 
sults of our subsequent procedures for the separate diagonals, and a 
steady-state for the whole circuit is attained only when one has ac- 
crued for each separate diagonal. 

We may define a set of intervals for every synapse $; , which we 
shall call the ranges of succeeding synapses s;.; from s;, and denote 
by RF; (si), by the following recursion: 


(1). If a; =.0. he Ss) eas oo). 
If Qa; < 0 wea (S;) — (cia ’ —o). 
(2). If KR; (s;) = (m,n), Rj (8i) = (oisja + MOjs , 


Cisjur + NA;,;), unless m or 7 or both are negative, in which 
case we replace m or n respectively in this expression by 0 . 


If the range of s;.; from s; is (m, ), the end-points m and n 
are denoted by L;(si) and U;(s:). If L;(si) = U;(s;) for some PO 
is said to be inaccessible from s;; otherwise accessible, Clearly, if the 
suffixes 7,7, h, 1 are in cyclic order, Ry; (8;) is a proper subset of the 
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range of s, from s;; it is also inaccessible from s; or any preceding 
Synapse. Physically, that the range of s; from s; is (m, ”), where 
say m < , means that by properly varying Yi , We can cause y; to as- 
sume any value within the interval (m, ); but we cannot from 8; 
cause y; to diminish beyond m, nor to exceed n. If S; 18 Inacces- 
sible from s; , then nothing happening at s; will have any effect upon 
the excitation at s; , and consequently, if a circuit contain any synapse 
inaccessible from anywhere, it will attain a steady state immediately, 
in the second n columns of E , and this steady state will be wholly in- 
dependent of the initial stimulation );. Incidentally, if at any time 
there are two zeroes between s; and s; , between which is an odd num- 
ber of inhibitory fibers, s; is inaccessible from s; (cf. Householder 
1941b, lemma 1). 

We shall now find it desirable to establish two lemmas of impor- 
tance, the first of which will be necessary for our more general theory, 
the second whereof is a special case of a result later to be established 
for general networks. 


LEMMA 1 

Let (m,n) be the range of a given synapse s; from s,. Suppose 
first that there are an even number of inhibitory fibers between s; and 
s;. Then we can find two numbers, 3 and 0, ,3) = 3, such that for 
ali values of y,; = 3, the consequent value of y; ism; ford, <yi Sh, 
the chain between s, and gs; is completely active, and y; is accordingly 
a linear function of y,; for yi > 01, y; isn. Second, suppose there is 
an odd number of inhibitory fibers between s,; and s;. Then we can 
choose-), and 9, , 3 =%,, so that fory; =%,y; =n; forme =yi <h, 
y; is a linear function of y; , and there are no zeroes between s; and 
s;;foryi >h,yi=m. 

For convenience in statement, we shall assume an even number 
of inhibitory fibers between s; and s;.;; the proof may be extended to 
the other case with inessential changes. It is clearly sufficient to show 
(1) that as we increase y; from 0 to + o, then, unless L; (3 eet, 
when 3, = — o, that intermediaté zero between s; and s;,; which is 
nearest to s;.; is removed at most once, to produce a state of complete 
activity; and (2), if we augment y; still further, so that a new zero 
is formed—otherwise U;(s;) = + 0, = + co—this zero cannot be 
removed by further increase in y;. (1) follows thus. At values of Yi 
sufficiently small so that if R;(s;:) is finite yj+i is equal to it, there will 
be zeroes at synapses between s; and sj,; , of which one, say s; is clos- 
est to s;.;. Now suppose that, as we increase y; , the zero at S; is re- 
moved for a least value of yi, SAY Yio- Now if there were a zero be- 
tween s; and s, say at s,,, which remained unchanged as y; assumed 
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the value yj. , 8; could, of course, not be affected. If all zeroes between 
s; and s, disappeared for y; = Yio, a state of complete activity would 
result, in accordance with (1). The only other possibility is that there 
be no zero between s; and s; for yi < Yio, but one is initiated, say at 
s,,, for yi = Yio. Now if there is to be no interval of complete activity 
at all, then, precisely when the zero at s,, is formed, or y,, = 0, and 
not at all before, s;, must cease to be a zero, so that y, > 0. Now yx is 
a linear function of y,, for all values of the latter greater than or 
equal to zero, since there is then no zero between them: consequently 
if, as is clearly possible, we select a small y;, satisfying 0 << y¥; < yx, 
y’ will be the value of this function for some argument y’,, which 
exceeds zero, by continuity and monotonicity; and since y’, > 0, the 
zero at s, is broken down for this value of y; also; and consequently 
for ay, > 0, which is contrary to hypothesis. Hence (1) holds. (2) 
follows immediately: if we have complete activity between s; and 3; 
for a given range of values of y; , then all y, for i < k < 7 vary linear- 
ly with y;; and it is only when one of these, say y;, diminishes (lin- 
early) to zero that the complete activity is interrupted. Now if com- 
plete activity were to be reéstablished for a higher value of y; , then 
y,, being the same monotone function of y;, would be exceeded by 
zero, which is not compatible with complete activity. By (1), the 
zero formed at s; could not be removed in any other way. 

It is clear, by the continuity of the value of yi; as a function of 
yi, that the quantities 3; , 3;, and the ranges are connected by the 
important relations 


L;(s;) — A iis; Bo + oi”, 
U;(8i) = Aijing A. + ing ™. 
For the pair 9), 3, in the case where s; = Sim = s;, as determined 


thus, we shall use the permanent notation %;, 0,;. In this case we 
have 


Ln(8i) =A Voi + ui, 
U,n(si) =ABi + Wi. 


LEMMA 2, 

I i a circuit C has been in complete activity for a period qn time- 
units in length, and § <1, then the excitation of s, at the time qn + & 
is given by 


1— As 
yilqn + 6) m5 +a Ae, (1) 
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are that at the synapse s,,;,j <n, for any time qn+j+é&,é&<1, 
y 
AQ 


Yirs (QM + I + E) =i + Ain mi a pee 2 


If C is in complete activity we may write a difference equation 
for excitation at t + n as a function of that at t; this is 


yi(té+n) =uwtAyi(t). 


One verifies that equation (1) is a solution of this by substitution; 
and equation (2) is an immediate consequence of equation (1). 

Armed with these results, we shall easily be able to determine the 
possible types of activity in the circuit CG. To this end we may divide 
the possible values of the activity y; at s; into three regions: first, the 
region I’, , comprising those which are less than #; , second the values 
Vo < Yi = Vix, making up I, and third, the y; > #1, which consti- 
tutes I’;. Suppose the value of y; at a given time belongs to I’,, i.e., 
Yi = % i. By lemma 1, the consequent value of the impulse returning 
to s; att + n will be A %; + w;. If this is also less than %; , the same 
value will recur at t + 2n; and we shall have a steady state of C de- 
termined by y; = A %; + wi; if, however, A J + ui 2 oi, the circuit 
will, by lemma 1, go into complete activity, and there will be no steady 
state with y; = Ai doi + wi. 

Suppose that a given y; belongs to I2, i.e., that Io, S yi < Hui. 
Then C will be in complete activity, and, by lemma (2), the course of 
activity will be described by 


1 At 


,; — Mi + At A; 
Yi BL T—A 


\ 


Consider the value of this expression as t becomes very large. If 
|A| > 1, then, unless both yw; and 4; vanish, in which case also C is 
not in complete activity, y; will increase indefinitely, so that a steady 
state with complete activity is not possible. If A = + 1, we arrive at 
a special case to be treated later. If |A| < 1, yi will approach an 
asymptotic value 


fae 
Ys aaa? 


and if this value for y; puts C into complete activity, it will determine 
a possible steady state of complete activity: i.e., this will be the case 
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if and only if 
Mi 
| —A 
The possibility of a steady state for a y; in I’; may be treated ana- 
logously-to that of J’,: we shall have, as necessary and sufficient con- 


dition for the existence of such a value y; = A Hi + uw that A dG + 
ui 2 b1;. Collecting these conditions, we shall have the 


Doi = 


< Pris 


THEOREM. 
A circuit C may have at most three possible steady states: 


A). One given by y,;=ADoi + Mi 
pal 
1—A 

C). One given byy,=Ahi tu. 


Necessary and sufficient conditions for the actual existence of each of 
these are given respectively by the requirements: 


A). a= (1 — A) %; . 


B). One given by yi = 


Mi pe 
1s So EG AVE e 
B). Ba <5 =i) |A| 


C). hi > CA) Or: 
By way of corollary from these conditions—which differ from 
those of (Householder, 1941b and c) in giving explicit expressions for 
the equilibrium points—we may derive the general results of (House- 


holder, 1941b and c) for the simple circuit in a very much easier way. 
Suppose first that 0 <A <1. Then (A), (B), and (C) become 


Ui 


NA Ot ee Ge. 

ADT oe 

B). 9: <7—GS¥u, A] <1. 
Mi 

CA ere: 


Clearly one, and only one of these conditions can be fulfilled by a giv- 
en circuit. We shall consequently have the 
COROLLARY. 


When 0 < A <1, the circuit C has a unique steady-state activity 
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pattern. This will also be the case when A < — 1 , except that here 
no steady state of complete activity is possible. 
Suppose now that A > 1. Here again, condition (B) is excluded, 


and there is no steady state of complete activity; the other possibil- 
ities may be expressed as 


In oma 


Cys 


which are not incompatible, so that we may well have two possible 
steady states, given by A 9; + wu; and A 3, + mw; respectively, both of 
which contain at least one zero. 

It may be desired to trace the course of activity in C from the 
initial value 1; to whatever steady state, if any, is finally reached. 
This may be done without difficulty by considering which region con- 
tains 4,;, the initial value of y;. If this region contain a possible 
steady-state value of y;, then, if it is 7, or I’; the steady state will 
obviously be attained immediately, while if it is ., it will be ap- 
proached asymptotically, in accordance with (1). If the region con- 
taining 4; do not contain a steady state then, if it be 7’, or I’; the value 
of y; will enter immediately and move through the adjacent regions 
until the first one containing a steady state is entered, whereupon 
that steady state is either attained immediately or asymptotically ac- 
cordingly as this final region is I’, or l;, or is I’,; while if the initial 
region of y; be I’, then, if only one of I’, , f, contain a possible steady 
state value, then that one will be entered, at a time determinable by 
solving (1) for ¢ with y; = %; or 3; accordingly as the final region be 
I, or I'3; but if both I’, and I; contain possible steady states, y; will 
be described by (1) until that value y; = %; or 3,; which occurs first 
befalls, whereupon it goes into the steady state of I’, or I’; respec- 
tively. 

It will be worth while to conclude our discussion with a consider- 
ation of the interesting and important case where A = + 1. We no- 
tice first that if A = + 1, then the conditions (A) and (C) of the 
theorem become 


Ae?) : Mi = 0 
Git) : [i = 0 ; 


at least one of which must hold, so that there always exists a steady 
state for C with at least one zero in this case. The interesting case 
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arises when we have complete activity, however, at least initially, so 
that we shall suppose #; < 4; = 01: . Equation (1) becomes in this case 
| Oe ee 

Yar eed (aT 4» 
for the values of ¢ specified in lemma 2, which is indeterminate; evalu- 
ating the quotient by limits we obtain 


yi=t wi (+1)? + (4 1)1. (3) 


If uw; ~ 0, y; becomes indefinitely large with ¢, so that a steady state 
of complete activity is not possible; and the case has no special inter- 
est. When »; = 0, however, we have 


Ye Crd) as (4) 


In the case 4; = + 1, this means that the circuit C has an infinitude 
of possible steady states, one for each value of 4; satisfying 3) = di 
< %1:, so that the parameters of C do not determine the steady state 
conditions even among a finite set of values. Another such case arises 
when 4; = — 1; here, however, we have no steady state at all, but 
instead continual oscillations in C of amplitude |/;|. It may be re- 
marked that this possibility might have been determined by the meth- 
ods of (Householder, 1941a, b, c) except for an oversight: the appli- 
cation of Cramér’s rule on page 68 of (Householder, 1941a) is invalid 
when A = 0, which is equivalent to supposing A = 1(4—=1-— A);on 
account of this the principal result of (3) is not in general valid. If 
the = sign be struck out of the statements of this result, however, and 
replaced by < , it becomes correct. A similar correction can generally 
be made at a few other places in (Householder, 1941a, b, c, 1942) 
where this.oversight enters. The two possible cases in respect of the 
equations (7) of (Householder, 1941a), namely inconsistency and in- 
determinacy, correspond respectively to the case where u; does and 
does not vanish for alli. In the first case, as we saw above, equation 
(2) gives an indefinitely great y; for sufficiently late t , and no steady 
state of complete activity is indeed possible; in the second we have the 
determination of the substantial one of an infinite number of possible 
steady states by the initial values 4;, as we should infer from the 
equations (7) of (Householder, 1941a) in this case. 

We may perhaps conclude the present discussion with two re- 
marks. First, with regard to the relationship of our present analysis 
to the purely static one of (Householder, 1941a, b, and ¢c) ; we have 
been able to obtain all ‘the results there found with a rather simpler 
procedure. We are, moreover, enabled to solve explicity a problem 
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treated less directly in (Householder, 1941a, b, and c). A theorem 
there enables us to calculate all the stimulus patterns consistent with 
one or two specified activity patterns; since there are only a finite 
number of activity patterns, we can always enumerate them until we 
find all the pairs consistent with any given stimulus pattern, which is 
what we really want. This method may be very laborious, however: 
there are 2” distinct AP’s, and 2”.(2” — 1) pairs to be tested in gen- 
eral. Theorem 1 above on the contrary enables us to determine at once 
explicitly how many steady states there are in any given circuit, and 
what values of y; determine them. 

Secondly, the above results may be extended (1) to the case 
where the initial stimulation between t = 0 and t = 1 is not a con- 
stant 2; , but a bounded function y;(t), by simply dividing the inter- 
nal (0, 1) into sufficiently small segments so that all points of a given 
segment, when taken as the J; above, permit the circuit to approach 
the same steady state in the same way; the predictions thus obtained 
for each segment are to be applied only throughout those later time- 
intervals which are congruent to it modulo n; (2), to cases where 
the initial variation in stimulation has lasted arbitrarily long before 
becoming steady, by waiting until the last synapse of (C receives a 
steady stimulation, counting that interval as the first, and applying 
(1). 
Desinently, I should like to express my appreciation of the coun- 
sel and criticisms of Dr. A. S. Householder, who suggested the prob- 
lem of this and later papers on the theory of neuron networks. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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DELAYED ADSORPTION AND DIFFUSION IN 
COLLOIDAL MEDIA 


HERMAN BRANSON 
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The behavior of the diffusion coefficient of a solute which can be ad- 
sorbed by a colloid only after the colloid has aggregated to a certain size 
is deduced on the basis of a few assumptions. Some relations of such a 
mechanism to cell reactions are indicated. 


A type of reaction plausible in cell behavior is one in which two 
or more substances may be released under a stimulus. This paper 
considers such a release with two substances one of which can aggre- 
gate into larger micelles and upon reaching an optimum size adsorbs 
the second substance, which we call the solute. After this initial ad- 
sorption we may consider either no further aggregation and adsorp- 
tion or further aggregation with adsorption as being the possibilities 
of greatest interest. For our purposes the significant point is what 
effect these alternative possibilities have upon the diffusion coefficient 
of the solute. In addition, we indicate how this mechanism may be 
used to interpret certain reactions in a cell, e.g. reactions which begin 
at a certain rate, proceed at that rate for some time, and then fall to 
a minimum. 

The aggregation of the colloid particles is assumed to take place 
in the following chain: 


Ny, + Ny 7 No, 


by ae We a) 


(1) 


Nm a i Ty 7 Nin , 


where 7,, is the size at which adsorption of the solute occurs. This 
chain expresses the assumption that a higher aggregate is formed 
from the next lower by the adjoining of a simple micelle, in short 
there is no aggregation of higher aggregates with each other. We 
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shall consider first the mechanism which allows no aggregation after 
adsorption. The differential equations for this chain are 


dn, 
Pe oe k, N,? SE ky Ty Ne =e Kim Ny Nm » 
dn, kin? 
——= —khn, nN. ; 
dt 9 2 2 
2) 
dne 
af = Ken Ny Near — tee Ny Ne, 
TNin 
‘dt ose Mma Ny Ny. 


On making the substitution n, d t = d x these equations are trans- 
formed into the linear forms 


dn, — m-1 
dx << ~ ki Nj r) 
dn, ey 
dz a 2 Ny key Ne, (3) 
dn, 
Aa = Kn-a Nn-s- 
The characteristic equation of this set is 
A+ ky ky Lego rss cate ine 0 
k 
ae ol A = ks 0 S ihe tt ee 0 
=0. (4) 
0 Pasi ie A ate k, © elites oe bee celts 
SPU at eet emia en a ee che tel Be Cat Oc Chruke Fix el a 


A simple relationship is found between the determinants of each or- 
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der. Consider for m = 8 we have 


Asics ke 
D,(4) = k, 5 
Seige Anetra hon 


observe that 


D;(A) = (A + kz) Do (A) oe 
and in general 
k, ka +++ ke 
D(a) = (A+ Be) De-a(A) + pee pee(5) 


Inasmuch as we shall not treat the general case, except to observe 
that the physical conclusions will probably not be much different for 
large values of m, we can state that in order to be physically mean- 
ingful the solutions of equation (3) 


m 
Ny, = Dd Cruz er (6) 
j=1 


satisfy the boundary conditions at all values of time 
3 k Nk; = No. 
1 


Since x and not time occurs in equation (6), we see that it is not nec- 
essary that the real part of the 4 be negative. From the defining equa- 
tion for x, we have x > %, %o finite, n, > 0, then t(7%) > oo; thus 
an infinite time is required for all the n, to disappear even though the 
Xo is finite. 

All the information needed for our purpose can be had from a 
detailed treatment of the set for m ='3. That is there will be an 
aggregate of three colloid particles built up before adsorption of the 
solute occurs. The problem immediately suggested is the relative be- 
havior of n, and n.. This behavior can be obtained from the integral 
curves. Using the notation of L. R. Ford (1933), 


ed Bs + ky Ne 
AL Ewes ay Bee (7) 
dn, ky Ny + Keg Ne : 


A= (ky + x)? -6hi ke. 


Since k, + k, # 0, the integral curves are not conics. 
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. = 2 
M,: Nz = [hy My Moe Ns — Ub Ny 


_ -— (4h) + v4 


SE petemmemrenG vier ee (8) 


the integral curves are shown in Figure 1. 


FIGURE 1 


These solutions are physically well-behaved in that all show that 
either n, disappears before n. or they disappear together. But the 
presence of n, after the disappearance of 7, , makes it impossible for 
all of , to be transformed into n;. Thus for the simple case of n,, 
N., N23, we can expect as a resultant state only 3; or a mixture of 7. 
and n;. The first will result if k, > k, and A 20. The k’s will be de- 
termined by the reaction, hence in a specific experiment we can ex- 
amine the products present after sufficient time which will decide 
whether n, and 7; or only n; remains. 

An alternative treatment of the integral curves which is imme- 
diately applicable to the general case has been suggested in a letter 
by Dr. A. 8. Householder. From the conditions, the integral curves 
in the n,, m, plane are confined to the region bounded by the coordi- 
nate axes and the line n, + 2n. = n,.. Since dn,/dx < 0 for any n, 
and 72 not both zero, the integral curve can cross the n. axis (n, = 0). 
While for n, = 0, unless n, = 0 simultaneously, the integral curve is 
directed into the region. The scalar product of the vector with the 
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normal is —3k.n, which is always negative for n, > 0. Initially 
NM, = 0 and n, =n, , thus the curve begins at the intersection of the 
line n, + 2n. =.) and the », axis; the integral curve is tangent to the 
line and directed upward from the n, axis. As n. builds up the inte- 
gral curve turns inward away from the line n, + 2n,=n,. 

The general set is treated in the same manner. Here the integral 
curves are confined to the region bounded by the m—1 coordinate 


m-1 
hyper-planes and the hyper-plane } k n, = n,.. The integral curve can 


cross the hyper-plane , — 0 for always dn,/dx < 0. But when ™ 
=0,kF-1, dy/dze=— hy. m, > 0 unless % is also 0; but if m. is 
0, dn,./dx, > 0, etc., until we meet an n, # 0. Hence if any 7, 
k #1, approaches zero the integral curve turns toward a region where 
the'n, is increasing. Finally the scalar product of the tangent and the 
outward normal is —(m—1)km-1 %m+ which completes the proof that 
for a finite x , but an infinite t, m, vanishes in the general case; and 
although some of the other 7’s may vanish simultaneously with 7,, 
they cannot vanish ahead of 7, . 
The solutions of the set satisfying the boundary conditions 


N, + 2n, + 8n3; = N, 


and at x = 0, n, = % where % is the initial number of simple colloid 
particles, are 


n,=—* | (a + ky)e* — (Ap + ken] 


A A> 
BFE H DS oe], 
kt (As — 2) (9) 
No 3(A1 aE k,) (he + 1) ee 
Le ae Zl Daas) 
AaCAy —s ho) 


where _ 
—(k, +h.) + VA 


daa = 5 


The relation between x and ¢ is ™% t= 2 


= [= 
euls 
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where x = 0 when t= 0. Before the integration can be performed 
the upper limit must be fixed: Calling this upper limit # and.recalling 
that we are concerned about the value of the time from the beginning 
of the reaction until the appearance of the first. particle capable of 
adsorbing a solute particle, we have on introducing n; = 1 into (9) 
and calling the coefficients of exponentials a and f respectively 


3/% —1+ae— pert (10) 

This equation cannot be solved explicitly for «. Making the substi- 

tutions ; 0) 
Yi=—aer+1—3/n, 


Us = B exe ¢ j \ 4 
we can solve graphically if we have numerical values for a, $, and 
my’. Taking 753-1000; e=4 k= 4k ithenyA = Myedg a es 


4. = — 6k and we find a = — 5.55, 6 = 3.75. Substituting these 
values we have x = 0.48/k 


= —_ ——_—— dx, 

1000 J Dio ce 

where 7 is the time required for the first n, to appear in the system. 
Integrating numerically by using Simpson’s rule, we find 


7 = 48.18/1000 k secs. » (11) 


Our analysis culminating in (11) can be summarized: when the 
aggregation takes place according to (3) with m.= 3, at the end of 
48.18/1000 & secs. the particles capable of adsorbing solute particles 
appear in the system. Before their appearance, the other colloid par- 
ticles would have ‘negligible effect upon the diffusion coefficient of the 
solute. After their appearance, however, they would act to reduce 
the value of D according to J. Reiner’s curve (1939). The effect is 
Shown in Figure 2a. If the aggregation continues with the number 
of adsorbed solute particles a function of the surface area as in an 


DIFFUSION COEFFICIENT 
OIFFUSION COEFFICIENT 


2 


TIME ; i TIME 


FIGURE 2 
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earlier paper (H. Branson, 1942), we would have the variation rep- 
resented in Figure 2b. 

In order for the mechanism to be applicable to cellular reactions, 
7 has to be in general small, except for some reactions where it may 
be of the order of a second. From physical considerations we see that 
k is a function of the diffusion coefficient and the radius of the agere- 
gating particles; M. V. Smoluchowski (1918) found for a colloid so- 
lution k = 427 D. Introducing values given in that paper from 
Zsigmondy’s experiments we have k ~ 10-2 and + & 10?° secs. Thus 
unless the diffusion coefficient within the cell is considerably larger 
than in solution, this reaction gives an inordinately long time for the 
beginning of the decrease in the diffusion coefficient. More plausible 
values of 7 can be obtained by considering 7, to be much larger in 
(10). Raising n, to 10’ causes practically no change in x, and taking 
r © 10-* cm then for D of the order of 10?, 7 will be of the order of a 
hundredth of a second. This is not an unreasonable value of D for 
aggregations where the binding energy is large. 

The author wishes to express his thanks to the Julius Rosenwald 
Fund for a grant which supported the initial steps in this problem. 
His thanks are again extended to Prof. N. Rashevsky and the Uni- 
versity of Chicago group of mathematical biophysicists, especially Dr. 
A. S. Householder, for discussions and suggestions. 
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AN EXPRESSION FOR THE RATE OF RETURN OF AN EGG 
AFTER ARTIFICIAL DEFORMATION 


H. D. LANDAHL 
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The forces which may be involved in the restoration of a deformed 
cell to its normal shape are considered. Estimates of the order of magni- 
tude of the forces suggest that the most important forces are those due 
to surface tension, membrane elasticity and viscosity. An approximate 
expression is then derived for the rate of return of an elongated or com- 
pressed egg. The former expression is compared with data on eggs of 
Arbacia by E. N. Harvey and H. Shapiro, and it is found to agree suf- 
ficiently well with the data. The initial surface force is too low compared 
with that given by K. S. Cole but various factors are discussed which 
could contribute to this discrepancy. The value for the net viscosity of 
the egg is about twice the value arrived at by L. V. Heilbrunn for the 
egg protoplasm. 


If a cell is deformed by an external force, various forces may 
tend to restore it to the normal shape. These may be grouped into 
surface forces and volume forces. The surface forces may be sub- 
divided into surface tension, elastic membrane tension and osmotic 
forces. The former is independent of area, while the second is prob- 
ably directly proportional to the change in area, at least over the 
initial range. The volume forces may be due to electrical or diffusion 
phenomena. The effect of electrical forces can probably be neglected 
(Rashevsky, 1938, Williamson, 1939). The same may be true for elec- 
trical surface forces. If the period for return to normal after defor- 
mation is short compared with the time constant, A, of diffusion pro- 
cesses (Landahl, 1939), we may consider the diffusion field as con- 
stant. For Arbacia eggs, A for oxygen is about four seconds, and that 
for other metabolites probably considerably greater. Then the effect 
of diffusion and osmotic forces will have little effect on the form of 
the curve of change of shape. As the eggs do not appreciably deform 
under the action of gravity, we shall neglect this volume force. Thus 
none of the volume forces will be considered. 

Opposing the surface forces are those due to the inertia and the 
viscosity of the material of the egg itself, the latter being the result- 
ant of the viscosity of the protoplasm including the effect of the gran- 
ules, internal structure, and membrane, as well as the viscosity of 
the surrounding medium. The effect of this latter viscosity is roughly 
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to add to the net effect of the cell viscosity. This may be seen in the 
expression obtained by J. A. Wheeler for the “relaxation” time in the 
case of a true surface tension (Harvey and Shapiro, 1941). Since the 
effect of the cell is then at least several times that of the surrounding 
medium, the assumption of simple additivity introduces,.a negligible 
error. The inertial terms are not more than a fraction of a per cent 
of the term due to surface forces. They may be estimated as follows. 
At any instant the force due to inertia will not be of a greater order 
of magnitude than if the half mass were accelerated at the rate at 
which the pole is accelerated. This product of the half mass with the 
acceleration gives a force which may be compared with the effect of 
surface tension for that particular degree of elongation. An estima- 
tion indicates that the former is at most a fraction of a per cent of 
the latter. An assumption of an elastic membrane force only increases 
the effect at greater elongations and it is only for these shapes that 
the inertial effect is at all appreciable. Hence we shall not consider 
the effect of inertial terms. 

Let 27, be the length and 2r, the width of a deformed cell whose 
diameter in the spherical shape is 27,. We shall assume that any de- 
formations can be approximated by an ellipsoid of revolution about 
the z-axis, along which 7, is measured. Let the surface force per cm. 
at any point be y’( x,y, 2) = y'(p, 2), where p? = 2? + y?. If Ki 
and K, are the principle curvatures at (p, 2) and y; and y’, are sur- 
face forces per cm. in the respective directions, then the normal com-_ 
ponent is y’,K, + y’.K,. The component of y’ which is due to a sur- 
face tension is independent at deformation, and can be treated with-" 
out too great difficulty. The component due to an elastic membrane 
tension is of a different order of complexity. For this reason, in or-' 
der to treat the problem, we shall introduce an approximation and for' 
the present only consider the deductions which then follow. We shall 
assume that the elasticity is the same in any direction, though, in 
general, it may not be isotropic. We shall introduce an average sur- 
face force y which is the average value of y, + y2 over the surface for 
a particular deformation. We shall further assume that is varies ap- 
proximately linearly with increase in area, or that = 


CAS {Ea k 
y=ro+ 3 ; ty Adda 


where S, should be taken as the surface of the sphere of no membrane. 
tension, 4 S is the change in the surface and € is. a constant.: With-:, 
out much error one can let S, = 4.2 7,2? and measure the area changes . 
from, this value. That the relation of: (1) is then approximately lin-: 
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‘ear in the case of Arbacia eggs is indicated by unpublished data fur- 
nished by Dr. K. S. Cole. The value of £ is in this case about 0.58 
dynes per cm. If one takes the value for the membrane tension .19 
‘dynes per cm. given by Harvey (1931) for a twenty-five per cent 
increase in surface area and uses the value .08 for the initial value 
(Cole, 1932) one obtains ¢ = 0.44 dynes per cm., a value which agrees 
rather well with that above. 

It may be pointed out in this Baer that the surface force 
‘measured by Cole (1932) by the compression method is a quantity_of 
the same type of average as y used here. Evidently when the shape 
is only infinitesimally different from a sphere the error in y is neg- 
ligible. But for large deformations the error introduced may become 
considerable. 

The relative rate of change of the half length, 7, , under the in- 
fluence of surface forces only, is given by Betti’s theorem applied to 
plastic flow by the seep 


1 adr, _ 
ta at ote) nV 


—— ff [22-— 47 Xe yy Y) Ids, (2) 


where V is the volume of the cell, n is the viscosity, X , Y and Z are 
the components of the surface forces along the x, y and z axes. The 
integration has been carried out by G. Young (1939) for the case in 
which the cell is a prolate ellipsoid of revolution. The result is 
[Young, 1939, equation (72) ] 

1 dr, y (1% —4'o) 

rm dts on n Civ (3) 
where 6 = 1 for 7, =7 2, 8 =.94 for 7; = 27, and 6 =.74 for 7, > > 7. 
We shall be interested in a range 7, — 7, to 7, © 27, so that we may 
set 6 = 1 without appreciable error. A consideration of the case of 
the oblate ellipsoid shows that equation (3) again holds, but 6 takes 
on slightly different values. Considering constant volume so that 


1412? =1 >’, (4) 


and setting 


pee te (5) 


Lo 


we may write equation (3) in the form 


53 
Sets e)9?— 1), (6) 
dt NY o 
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In order to integrate equation (6), it is first necessary to obtain 
the relation between y and «. This follows directly when we find the 
relation between AS/S, and «. Substituting the value of the area of 
a prolate ellipsoid of revolution into (S — S,)/S. and using equations 
(4) and (5) we have an expression for 4S/S, in terms of «. Evaluat- 
ing the-square root of 4S/S, numerically for values of ¢ and plotting 
the result one finds an approximately straight line relationship in the 
range « = 0 to e = 0.6. Deviation becomes marked beyond this, the 
function being concave downwards. For the oblate ellipsoid, the curve 


becomes markedly concave upwards beyond about « = —0.4. In the 
range of « from —0.3 to 0:6 we may write 
= 04 6(1—0.86 ++), (7) 
0 


where the error is about twenty per cent at the extremes. The range 
in de/dt which we wish to consider is found experimentally to vary 
by a factor of the order of a hundred in the range of from 0.1 to 0.6. 
Thus a considerable error in de/dt at any point can be tolerated. 

Introducing equation (7) into (1) and the resulting equation in- 
to (6) we obtain, 


de__ 
di 
where we neglect the factor (1 + ¢/4 — «?/24 + ---) from the ex- 
pansion of the brackets in equation (6), as well as the factor (1 — 
0.8 « +---) of equation (7). The error in neglecting the first factor is 
not more than fifteen per cent, while that of the second factor may 
be about fifty per cent at « = 0.7, but in the opposite direction. It is 
equivalent to requiring that 0.4 ¢ is larger than y.(6/e — 1)/24. 
Since we shall not consider values of « less than about 0.05, this 
amounts to requiring that y, should be less than about 0.05. 
Integrating equation (8) and setting « = «, for t= 0, we obtain 


63 7 To Eo" (yo + .4€ €?) 
t= 05 
Yo 7 (yo + AC Eo") 


The above expression might be expected to apply approximately 
for the recovery of a cell which had been compressed between two 
parallel planes and then released, as well as to the case in which a 
cell was elongated as a figure of revolution and then released. The 
necessity to assume an averaged y in the integration presents itself 
as a difficulty. The error involved may not in general be the same for 
the prolate and oblate ellipsoids. Thus, to that extent, agreement be- 


a ea (yo + 0.4¢ e)e, (8) 
"To 


(9) 
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tween corresponding parameters obtained from the two procedures 
would tend to justify the approximate treatment used. One would 
allow a slight overestimation of perhaps twenty per cent in one case 
and a corresponding underestimation in the other due to the fact that 
the errors in equation (8) are in the opposite direction. 

It may at times be convenient to use the ratio of the length to 
the width as a measure. If this is designated by R, one may substi- 
tute for « the value « = R*/? — 1 in equation (9) and obtain an ex- 
pression for the ratio of length to width as a function of time. 

We shall now compare the relationship described by equation (9) 
with the experimental results obtained by E. N. Harvey and H. Sha- 
piro (1941) from Arbacia eggs returning to normal after having been 
passed through sufficiently small capillaries to elongate them. From 
the length of the cell, 7, , at any instant, the value of ¢ can be ob- 
tained by equation (5), so that an empirical ¢(t) relation is obtained. 
Comparing equation (9) with this relationship, one can determine 
the parameters. The «(¢) relation is thus determined, and because of 
equation (4) and the condition for constant volume, 7,(t) and 7. (€) 
are determined. 

In comparing equation (9) with the empirical relation, it is 
found that the value of y, is required to be less than about 0.01, (and 
greater than about —.002) or the empirical relation does not even ap- 
proximately agree with the theoretical. Within this range the ratio of 
the parameters 7/¥y, to 6/y, is nearly constant and approximately the 
same as that obtained when y, is taken equal to zero, for which 
n/€ = 0.36 cm. sec." Hence, we shall use e(t) as if y, = 0, though 
only requiring that it should not be larger than about 0.01. This re- 
striction is sufficient to satisfy the previous condition that y be less 
than 0.05. Instead of equation (9) we then obtain the equation 


U/] Lh ee 
= |—_—___.. 10 
JS ey 


Introducing this expression for e« into (4) and (5) we obtain 7, (¢) 
and 7,(t). Since it is the length, 27, , and width 2r, , which are meas- 


ured, we write 
Qr, = 2ro[1 + {e02/ (1 + .64 ¢ en?t/n 7) }*) , (11) 
Qr, = 2ro[1 + {e02/ (1 + 646 &? t/y 10) }*) 3. (12) 


When 7/f = 0.36, & = 0.72 and 27, = 76 X 10+ cm., we obtain the 
curves shown in Figure 1. Using the value ¢ = 0.58 given above, we 
then obtain 7 = 0.21 poise. To obtain the net viscosity of the egg, 
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AND H. SHAPIRO 
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one should subtract approximately .01 because of the surrounding 
medium and perhaps allow for overestimation as mentioned above, 
giving 7- == 0.16 poise. The value 7- should be the total effective vis- 
cosity of the egg as a whole, if one could rely upon the approximations 
made. According to Heilbrunn, (1928), the viscosity of the fluid por- 
tion of the Arbacia egg is about .03 poise, while the viscosity uncor- 
rected for granules is about 0.07 poise. The presence of larger struc- 
tures may contribute appreciably to the total effective viscosity. Just 
to what extent this may be, would require a more detailed analysis. 
There is also the possible effect of the various membrane layers 
which affect the net viscosity. If we could assume that the value of 
ne = 0.16 is correct, then we should conclude that the contribution to 
ye by structures and membranes would be an amount approximately 
equal to that due to the protoplasm with granules. This would not 
imply that the amount contributed by the membrane gives its vis- 
cosity, though it may be possible to estimate it. In any case the value 
of 7. obtained is sufficiently plausible as to the order of magnitude. 
Using the value of ¢ from the data by E. N. Harvey given above, we 
obtain 7. = 0.11 poise. 

If the external viscosity instead of being that of wiltet were 
increased by a factor of 15, the above estimates would predict that 
the constant of the curve, 7/¢ , would be about doubled. Conversely 
one could thus determine the effective viscosity of the egg approxi- 
mately as that value of the external viscosity which ht sonbies the 
constant 7/¢. 

~ Thus far we have only required that», be of the fieibe or tests 
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than about 0.01 dynes per cm. The value for summer eggs, as meas- 
ured by K. S. Cole (1932), is 0.08 dynes per cm. Only one-half this 
value is obtained for late summer egg’s. (Cole and Michaelis, 1932). 
In the experiments under discussion, late season eggs were used. 
This is still probably too large a discrepancy to overlook, so we seek 
some explanation. Until the more exact solution to the problem, as 
mentioned above, is given, one cannot exclude the possibility that 
there would be no disagreement if the complete solution were known. 


Since equations (11) and (12) represent the data, we can say 
that the rate of change of ¢ is approximately proportional to «*. If a 
true surface tension is assumed, we see from equation (6) that the 
rate is approximately proportional to «. The approximation made 
initially, of an average effective y, probably has the effect of increas- 
ing the restoring force for large elongations. Also as the shapes di- 
verge from the elliptical shape in the direction of a truncated cylin- 
der, the effect of elongation is perhaps appreciably accentuated. That 
the shape is not quite elliptical is seen from Figure 1, where the meas- 
ured widths are somewhat lower than that expected, although the 
lengths are about correct. 

The naked surface of Arbacia eggs is apparently liquid, as indi- 
cated by experiments on coalescence (Chambers and Kopac, 1937). 
The outer layers are more solid in nature (Chambers, 1940; Kopac, 
1940). However, it is likely that a certain amount of plastic flow may 
take place in the surface. This results in a hysteresis phenomena as 
is observed in protein membranes (Harvey, 1936). This property 
would tend to reduce y, to some extent. But it is likely that the meas- 
ured value used above includes such an effect. Another factor not 
taken into account is the possible effect of volume elasticity (Pfeiffer, 
1937). However, E. Howard (1931) finds no evidence of plastic flow 
in Arbacia protoplasm over a considerable range of temperature. 
These observations did not include the firmer layers in the region of 
the surface. In general, if any part of the protoplasm has an effective 
finite yield point, the stable shape after a deformation beyond this 
point would deviate from spherical depending on the value of the 
elastic coefficient and the yield point. There is insufficient data at 
present to make an estimate of such a possible effect. 

As was mentioned previously, the time constants for the diftu- 
sion processes are probably well below the relaxation period. How- 
ever, the initial disturbance in the diffusion fields, due to the deforma- 
tion of the cell as well as to the presence of the deforming mechanism, 
might result in a force which tends to keep the cell slightly deformed. 
Thus the cell might approach a shape which is not a sphere (House- 
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holder, 1941). The latter shape would reduce to a sphere at a rate 
depending upon the time constants, A, of the metabolites. 

In view of the above factors, one cannot rely on the estimated 
value of y- as given above, but its order of magnitude may well be 
correct. At least it is not necessary to assume such a large value as 
is indicated by the analysis on the basis of a surface tension (Harvey 
and Shapiro, 1941). 

The data on Chaetopterus eggs given by E. N. Harvey and H. 
Shapiro do not cover a sufficient range to determine a permissible 
range for y.. However, if we neglect its effect and use equation (10), 
together with the relation « = R*/* — 1, the data in the range given 
lie close to the theoretical curve for which 7/¢ is about fifty. The sur- 
face forces for Chaetopterus eggs are quite comparable as to order of 
magnitude with Arbacia (Harvey, 1931). If ¢ should have the same 
value in each case, the viscosity in this case would be about one hun- 
dred and fifty times that for Arbacia eggs. 

Summary. An expression for the length and width of a cell, 
which has been given an initial deformation, as a function of time is 
derived by an approximate method on the assumption that the sur- 
face of the cell may be elastic. The expression is found to agree with 
data from Arbacia eggs but the value of the initial surface force is 
not more than about one-third of that obtained by fairly direct meas- 
urements. Some of the factors that may be responsible for the dis- 
crepancy in the initial surface force are discussed. The value of the 
viscosity thus estimated is of the order of that of the fluid part of 
the egg as determined experimentally. 

The author is indebted to Dr. N. Rashevsky and Dr. A. S. House- 
holder for reading and discussing the manuscript. 
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LETTER TO THE EDITOR 


AN EXPLANATION OF SPASTIC WALK DUE TO 
ENCEPHALITIS LETHARGICA 


In 1925 the writer suffered from an attack of ENCEPHALITIS 
LETHARGICA; was paralyzed on the left side; apparently made a 
complete and rapid recovery. In 1935, after ten years of freedom 
from after-effects, sclerosis apparently began—a slight spasticity in 
walking was noted developing on the left side. The flexor muscles’ 
reflex time gradually became less than the tensor muscles’ reflex 
time, causing a “kick-up” in walking. This has grown to be quite 
noticeable and does not disappear with the Belladonna treatment for 
after-effects. 

The kick-up can be explained with a simple assumption: 


(Atz) WwW im (Atz)r 
(Atr)w  (Atr)r 


=C, aconstant = 1, for normality. 


Here 
(At;) w= reflex time for extensor muscles at a walk. 
(Atz)e — 73 66 66 (73 66 at a run. 
(Atr) w = reflex time for flexor muscles at a walk. 
(Atr)r = 6é 6é 6c 6é 6é at a run. 
Clearly 
(Atz)w > > (Aty)e 
and 


(Atr)w > > (Ate). 
Then from the first assumption: 


and 
(Atz R— C(Atr)e A 


Let us define Aw as the difference in reflex times for the extensor 
and flexor muscles at a walk; Az as the difference in reflex times for 
the extensor and flexor muscles at a run. 

Then 
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Aw = (Atz)w — (Atr)w = (C—1) X (Ate) w 
and 
Ap = (Atz)r — (Atr)e = (C—1) X (Abr) e. 


But Ay is a measure of the kick-up at a walk, and Arg is a meas- 
ure of kick-up at a run. By inspection it is obvious that Jy > > Ar. 
This is precisely what I observe, for although I have a pronounced 
spastic walk, when I run the kick-up vanishes. This result, shown by 
my simple analysis, is a source of surprise to all who observe my 
running, yet it affirms the correctness of my simple first assumption. 

Pasadena, Calif. 
March 28th, 1942 
FREDERICK R. HIRSH, JR. 
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